Effective Kondo model for a trimer on a metallic surface 
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I consider a Hubbard- Anderson model which describes localized orbitals in three different atoms 
hybridized both among themselves and with a continuum of extended states (ES). Using a generalized 
Schrieffer-Wolf transformation, I derive an effective Kondo model for the interaction between the 
doublet ground state of the isolated trimer and the ES. For an isoceles trimer with distances a,l,l 
between the atoms, the Kondo temperature is very small for I < a and has a maximum for finite 
/ > a when a is small. The results agree with experiments for a Cr trimer on Au(lll). 



The study of many-body phenomena in nanoscale 
systems has attracted much attention in recent years. 
Progress in nanotechnology has made it possible to con- 
struct nanodevices in which the Kondo physics is clearly 
displayed [1-3]. Scanning tunneling spectroscopy has 
made it possible to probe the local density of states near 
a single impurity and Fano antiresonances have been ob- 
served for several magnetic impurities on metal surfaces 
[4-7]. These antiresonances observed in the differential 
conductance, reflect a dip in the spectral density of con- 
duction states caused by the Kondo effect [8]. Further- 
more, corrals built on the (111) surface of noble metals or 
Cu have been used to project the spectral features of the 
Fano-Kondo antiresonance (FKA) to remote places [5,8]. 
The observed Fano line shapes for one magnetic impurity 
on these surfaces have been reproduced by many-body 
calculations [8-10]. 

The situation is however different for Cr trimers on 
the Au(lll) surface [11]. Experimentally, no FKA is ob- 
served for a linear trimer, or a compact equilateral one 
in which the distance between any two Cr adatoms is 
believed to be the same as the Au-Au nearest-neighbor 
one (a = 2.88 A). For a single Cr adatom, no FKA is ob- 
served either. This implies that in these cases, if a Kondo 
effect originating from localized moments exists, its char- 
acteristic temperature Tk is less than the experimental 
temperature (7 K). On the other hand, surprisingly, for 
an isosceles trimer with one side of length a and two sides 
of length ~ 5.2 A, the spectrum does display a FKA with 
a half width Tk — 50 ± 10 K. This spectrum is similar to 
the one observed for one Co impurity on Cu(lll) [5,7]. 

The origin of this puzzling dependence of Tk of the Cr 
trimer with geometry is still unclear, in spite of recent 
theoretical efforts [12-14]. 

In this Letter, I start from a Hubbard- Anderson Hamil- 
tonian that describes three Cr adatoms interacting with 
a continuum of extended states (ES). Using a generalized 
Schrieffer- Wolff transformation (SWT) [18], the model is 
mapped into the usual one-impurity Kondo model. The 
resulting effective Kondo interaction averaged over the 
Fermi energy J depends on the geometry of the trimer. 
The ensuing Kondo temperature is consistent with the 
above experimental observations [11]. 

The form of the Hamiltonian is suggested by experi- 
mental evidence and previous theoretical work. The line 



shape observed for the isosceles Cr trimer with Tk ~ 50 
K is rather symmetric, while calculations using Wil- 
son's renormalization group for the degenerate Anderson 
model obtain a strongly asymmetric Kondo resonance 
[15]. This suggests that there is only one relevant 3d 
orbital. The surface potential breaks the degeneracy of 
the d orbitals into two orbital doublets with angular mo- 
mentum projection m = ±1 and m = ±2 and one orbital 
singlet with m = (the 3d 32 2_ r 2 orbital). Therefore the 
latter, singly occupied, is the candidate for the forma- 
tion of the Kondo state. A similar conclusion has been 
reached for Co atoms on Cu(lll) [8-10]. The absence of a 
Kondo effect for Cr [11] or Co dimers on Au(lll) at inter- 
atomic distances less than 6 A [16] suggests a strong in- 
teratomic hopping leading to a spin singlet ground state. 
Due to the symmetry breaking of the surface potential, 
the 3d 3 £2_ r 2 orbital of the adatom hybridizes with the 4s 
orbital, and then one expects a decay r _1 with distance 
r of the hopping between hybrid localized orbitals [17]. 
Finally the fact that Tk < 7 K for a single Cr adatom 
points to a small hybridization between localized and ES, 
justifying the use of the SWT. 

The Hamiltonian can be written as 

H = Hci + Hau + Hv i 
H C r = E d ^2 4(7 ^ + U 4t^»t4i^J- 

i<j 

- %(4<7 d i<7+ H - c -). 

k(T 

H v = J^e^'d^Cka + H.c). (1) 

Hcr is a Hubbard model that describes the isolated 
Cr trimer, with hopping that depends on the distance 
Tij = |Rj — Rj| between Cr adatoms as Uj = toa/rij. 
if Au describes a band of ES. I will discuss alternatively 
either bulk or Shockley surface states of the Au(lll) sur- 
face. Hy is the hybridization between localized and ES, 
assuming that the wave functions of the latter are plane 
waves [8]. 

If one eliminates Hy by means of a SWT neglecting 
Uj and uses another canonical transformation to map 
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the half filled Hqi into a Heisenberg model, the resulting 
effective Hamiltonian becomes analogous to previously 
used ones [12,13]. As shown below, this approach misses 
some coherence effects that reduce Tk strongly for the 
equilateral trimcr. 

I solve exactly Hq v retaining only the ground state 
doublet in the subspace of three particles 3cr) , and then 
I perform a generalized SWT to eliminate Hy . The va- 
lidity of the approach is discussed below. The result- 
ing effective Hamiltonian has the form of a one-impurity 
Kondo model 
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where S describes the spin operators of the doublet |3<r) 
and a acr i is a vector with the three Pauli matrices. The 
effective exchange can be written in the form 



Jkk' - %V V 
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where 
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\nsa) (\nta)) denote the six singlet eigenstates (three 
triplet eigenstates with maximum spin projection) of Hci 
with n =2, 4 particles, E v is the energy of the state \v), 
and E 3 is the energy of the ground state doublet |3ct}. 

The wave vector dependence of the effective exchange 
interaction J kk ' is rather involved. If only one virtual ex- 
cited state dominates in Eq. (3), approximating e k ~ e F , 
where is the Fermi energy (taken as zero), the effec- 
tive exchange has the functional form J kk ' = /(k)/(k')C, 
where C is constant. In this case, it is easy to see that 
if a Yosida trial wave function [19] is used to calculate 
the ground state energy, the result is the same if J kk ' is 
replaced by the angular average of J kk . Guided by this 
result I calculate this average evaluated at the Fermi en- 
ergy in the general case (calling it simply J) to discuss 
the dependence of Tk on geometry. It will become ap- 
parent that the qualitative features are not affected by 
this approximation. Except for some detail discussed be- 
low, I assume constant Vk- Then the expression for J is 
the same as Eq. (3) with the replacements 
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FIG. 1. Average exchange interaction J as a function of 
l/l for an isoceles trimer with two sides of length I and one 
side of length a. Full (dashed) line correpond to surface (bulk 
with ak,F = 1.6) ES. Other parameters are E d /to — —2.5, 
U/to = 5. 
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and e k - e F , with A? v = (3 | |dt»> % = ("l4«rl 3 T>. 
and Pij is the angular average of cos (k • (R, — Rj)) with 
e k on the Fermi surface. For an isotropic two-dimensional 
(2D) band [3ij = Jo(kprij), where Jq{x) is the zeroth or- 
der Bessel function and Ufo, — 0.6, The corresponding 3D 
result is /3,j = sa\(kprij) / '(kprij), with kpa ~ 3.5. How- 
ever, free electrons with e k = 6f and small wave vector 
parallel to the surface (implying small k- (Rj — Rj)) have 
a larger wave vector perpendicular to the surface. Nat- 
urally, one expects a larger extension of these states out 
of the surface and therefore a larger hybridization with 
the adatoms [8,10]. To simulate this effect I take kpa as 
a parameter smaller than or equal to 3.5 for the calcula- 
tions in which bulk states are assumed to dominate the 
hybridization with the localized states. 

In Fig. 1 1 show the average exchange J for an isosceles 
trimer with r 2 3 = a and m = r\ 3 = I, as Z is varied from 
+00 to a. I take £23 = as the unit of energy. Ed and U 
were chosen to correspond to an electron-hole-symmetric 
situation for only one adatom. Taking surface states as 
the ES, and disregarding small oscillations coming from 
Jq(x) at large argument, there is a moderate increase in 
J as the Cr atom 1 approaches the dimer formed by Cr 
atoms 2 and 3, with a maximum near I = 5a. For smaller 
I, J decreases monotonically until the equilateral shape 
is reached. The situation is qualitatively similar if the 
localized states are hybridized with bulk ES. The main 
difference is the larger value of k F - As kp is increased, 
the maximum is shifted to smaller values of I and larger 
values of J. For kpa — 3.5, the average exchange always 
increases with decreasing I, except for the small initial 
oscillations, while for k F a = 0, J is a monotonically de- 
creasing function [14]. 

To estimate the Kondo temperature for kpa = 1.6 
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(dashed line), I take a density of states per spin and site 
p = 0.15/eV, t = V = 0.5 eV, a band width D = 10 5 K 
[20], and use the expression [21] 

Tk = D^J~p~J exp ^^J^ • (6) 

This gives Tk = 8.3 K for the single adatom (I — ► oo), 
Tk = 9.7 K for the equilateral trimer, and a maximum 
Tk = 45 K for I = 2.3a. This behavior is roughly consis- 
tent with the experimental observations [11]. While Tk 
for the monomer and equilateral triangle are larger than 
the experimental bound Tk < 7, it is possible that taking 
into account more precisely the wave vector dependence 
of Jkk' improves the agreement. A scenario leading to a 
more drastic change in Tk is discussed below. 

For I — > oo, the average exchange coincides with the 
well known result for a single adatom as expected [18]: 
J = 2V 2 U/[(E d + U){-E d )}. The variation of J as I 
decreases is the result of mainly two competing effects: 
for (3i2 = 0i3 — 0, where there are no coherence effects 
between the atoms, the relevant matrix elements [Eqs. 
(5)] increase, while for I = a and fiij = 1, they cancel by 
symmetry. Specifically, if /3i2 = P13 is not too close to 
1, the largest contribution to J comes from the virtual 
state \2sg), the ground state of Hci in the subspace of two 
particles, which is a singlet even under the mirror plane 
of symmetry of the isosceles triangle m [22]. For I —>■ 00, 
only the ^4J„, B? v with j = 1 arc different from zero. 
As I decreases, the sums in Eqs. (5) increase due to the 
increasing contribution of A a - V , BJ V for j = 2, 3. Also, the 
energy level of adatom 1 is pushed up towards the Fermi 
energy due to its mixture with the bonding combination 
of adatoms 2 and 3 (the denominator E 2 sg — E3 + ep 
decreases). Both effects tend to increase J. 

However, there is a competing effect due to symme- 
try, as the compact equilateral shape is approached when 
kpl becomes small. For / = a, the Hamiltonian is invari- 
ant under the operations of the point group C% v . For 
kp — > — > 1), the operator entering Eqs. (4) be- 
comes J^i^lcn which is invariant under the point group 
operations. The state \2sg) is also invariant, while for 
I — > a, the ground state for three particles becomes four- 
fold degenerate with two spin doublets belonging to the 
two-dimensional representation E of Cz v . The matrix 
element of d iiy between any of these four states and 
\2sg) vanishes by symmetry [23]. Therefore, the contribu- 
tion of virtual states containing \2sg), which is the most 
important to J when /3y ~ 0, vanishes when fcj — > 1. 
This explains the decrease of J in Fig. 1 as I — > a [14]. 

The ground state for three particles |3er) is even un- 
der m for I > a and odd for I < a. As a consequence 
of this jump in the symmetry of the ground state, J be- 
comes much smaller for an isosceles triangle with one side 
longer than the other two, including the linear trimer as 
a limit. In particular, the contribution of |2s<jr) is negli- 
gible for I < a, since all matrix elements between even 
and odd states are very small. This agrees nicely with 
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FIG. 2. Same as Fig. 1 for surface ES and E d /t = -1.2. 
Dotted line is the difference cf + E2 sg — E^g. 

the fact that no FKA is observed for the linear trimer 
[11]. These symmetry arguments are independent of the 
approximation made in the averaging over wave vectors. 

The SWT is valid when the matrix elements connecting 
the ground states with excited states (Vk^4y' T , V^B^) arc 
smaller than the corresponding denominators. In addi- 
tion, since I retained only the ground state of in the 
three-particle sector, the effective matrix elements con- 
necting different states in this subspace should be smaller 
than the corresponding energy difference. For reasonable 
values of U, the eight states with one particle per site are 
well separated from those of higher energy. These eight 
states are degenerate if the separation between any two 
adatoms is large. However, in the situations of interest 
at least two adatoms interact, and then the spin quadru- 
plet is well separated from the two spin doublets and can 
safely be neglected. Precisely at the equilateral shape, 
the two spin doublets are degenerate and our approach 
is invalid. However a very small distortion towards an 
isosceles shape is enough to recover the validity of the 
SWT. This is because the distortion splits the spin dou- 
blets even and odd under m and the matrix elements 
between any even and any odd state are very small (they 
vanish for kpa — > 0). For the parameters used above 
to estimate Tk, a 1% change in the length of one side 
suffices. 

As mentioned before, the difference E 2sg — E^ + ep de- 
creases as I decreases from 00 to a, but is always positive 
and significant for the parameters of Fig. 1. However, 
for E d nearer to ep (for example ep — E d < t if U = 0) 
there is a crossing point and the ground state of He- 
rn, the equilateral trimer is the two-particle singlet \2sg). 
As a consequence of the decrease in the corresponding 
denominator in J, there is a dramatic increase in Tk 
for decreasing I near the crossing point, until the system 
enters an intermediate valence regime (where the SWT 
ceases to be valid) and for further decrease in I the sin- 
glet dominates the physics and there is no Kondo effect. 
The dependence of J on ! in a situation like this, taking 
surface states as the ES is illustrated in Fig. 2. The 
crossing point is at a/l — 0.576. Taking pV 2 /t n = 0.04, 
the resulting Kondo temperature given by Eq. (6) is 
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FIG. 3. Same as Fig. 2 for bulk ES, E d /t = -1.5 and two 
values of fcf- 

Tk = 0.33 K for the single Cr adatom, has a relative 
maximum Tk = 14 K at a/l = 0.2, and after a mini- 
mum near 7K, reaches 44 K at a/l = 0.47 and 174 K 
at a/l — 0.50. However, at this point the quantitative 
validity of the approach becomes questionable. 

Fig. 3 shows a similar situation when bulk ES states 
dominate the hybridization with the localized ones. The 
crossing point is at a/l = 0.75. The behavior of J is 
more monotonic than in the previous case. Taking again 
pV 2 /t = 0.04, one obtains T K = 0.05 K for I -> oo. For 
small afcp = 1.1, Tk reaches 51 K at a/l = 0.47 and 106 
K at a/l = 0.55. For ak F = 3.5, T K = 50 K at a/l = 0.52 
and 169 K at a/l = 0.55, where the SWT is at the limit 
of its validity. Beyond the crossing point (a/l > 0.75), 
no Kondo effect takes place. 

In summary, I have constructed an effective Kondo 
model for a Cr trimer on a Au(lll) surface. For an 
isosceles trimer with one short side and two longer sides 
of length I and reasonable parameters, the model pro- 
vides two scenarios in which Tk first increases with de- 
creasing I and then decreases as the equilateral shape 
is approached. In the first one, this decrease is due to 
symmetry effects and in the second one (which requires 
smaller ef — Ed), the ground state of the isolated equi- 
lateral trimer is a singlet and no Kondo compensation of 
a localized spin takes place. For an isosceles trimer with 
one side longer than the other two, or a linear trimer, 
Tk is very small due to the odd mirror symmetry of the 
ground state. While these results are qualitatively con- 
sistent with experiment [11], a quantitative agreement 
with all the different experimental observations seems to 
favor the second scenario, which is more robust under 
changes in the parameters, and that bulk states domi- 
nate the hybridization with the localized states of the 
Cr adatoms. Since in the second scenario the electronic 
occupation of each Cr atom in the compact equilateral 
trimer decreases by ~ 1/3, it is in principle possible that 
X-ray experiments [24] which are sensitive to the Cr va- 
lence can distinguish between the two scenarios. 
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